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Abstract 



o 

■ We study the evolution of quiescent galaxies at 0.5 < z < 2.5 as a function of stellar mass, using very deep 
NIR imaging data taken with the Multi-Object Infrared Camera and Spectrograph on the Subaru Telescope 
in the GOODS-North region. The deep NIR data allow us to construct a stellar mass-limited sample of 

\mS ' quiescent galaxies down to ~ 10 10 M even at z ~ 2 for the first time. We selected quiescent galaxies with 

Q^) , age/r > 6 by performing SED fitting of the multi broad-band photometry from the U to Spitzer 5.8^m 

bands with the population synthesis model of Bruzual & Chariot (2003) where exponentially decaying star 
formation histories are assumed. The number density of quiescent galaxies increases by a factor of ~ 3 
from 1.0 < z < 1.5 to 0.5 < z < 1.0, and by a factor of ~ 10 from 1.5 < z < 2.5 to 0.5 < z < 1.0, while that of 
star-forming galaxies with age/r < 4 increases only by factors of ~ 2 and ~ 3 in the same redshift ranges, 
t— » . At 0.5 < z < 2.5, the low-mass slope of the stellar mass function of quiescent galaxies is a ~ - 0.6, which 

is significantly flatter than those of star-forming galaxies (a ~ -1.3 - -1.5). As a result, the fraction of 
quiescent galaxies in the overall galaxy population increases with stellar mass in the redshift range. The 
fraction of quiescent galaxies at 10 n -10 115 M increases from ~20-30% at z ~ 2 to ^40-60% at z ~ 0.75, 
while that at 10 10 -10 10 ' 5 M Q increases from < 5% to ~15% in the same redshift range. These results could 
suggest that the quenching of star formation had been more effective in more massive galaxies at 1 < z ^ 2. 
Such a mass-dependent quenching could explain the rapid increase of the number density of ~ M* galaxies 

■ relative to lower-mass galaxies at z > 1-1.5. 

Key words: galaxies: evolution — galaxies: formation — galaxies :high-redshift 

1. Introduction quiescent population by quenching of star formation (Bell 

et al. 2004; Bell et al. 2007; Faber et al. 2007). Since re- 
. ' Determining how stars have been formed in galaxies cent studies suggest that the star formation rates in star- 
is crucial for understanding galaxy formation and evolu- forming galaxies seem to be simply almost proportional to 
*_h ' tion. In the present universe, it is known that galaxies can their stellar mass at least at z < 2 (e.g., Elbaz et al. 2007; 
5~ , be well separated into two populations, namely, passive- Daddi et al. 2007; Pannella et al. 2009; Kajisawa et al. 
evolving galaxies with red colors and star-forming galaxies 2010), how the quenching of star formation has occurred 
with blue colors (e.g., Kauffmann et al. 2003; Baldry et in galaxies becomes a key issue to understand the various 
al. 2004; Brinchmann et al. 2004). Bimodal color distri- star formation histories of galaxies. 

bution which consists of these two populations has also On the other hand, it is also known from studies of the 

been observed up to z ~ 1-2 (e.g., Bell et al. 2004; Wcincr 'fossil record' of present-day galaxies that the star for- 

ct al. 2005; Franzetti et al. 2007; Cirasuolo et al. 2007; mation histories of galaxies seem to depend strongly on 

Cassata et al. 2008; Williams et al. 2009; Brammer et al. their stellar mass (e.g., Heavens et al. 2004; Jimenez et 

2009). At z < 1, several studies have found that the stellar al. 2005). In the local universe, massive galaxies tend to 

mass density of passively-evolving galaxies increases by a have redder color and older stellar population, while most 

factor of ~ 2 from z ~ 1 to z ~ 0, while that of star-forming low-mass galaxies are young and actively star-forming 

galaxies evolves only very mildly (e.g., Borch et al. 2006; (Kauffmann et al. 2003; Brinchmann et al. 2004). Since 

Faber et al. 2007; Pozzetti et al. 2010; Ilbert et al. 2010). the pioneering work by Cowie et al. (1996), many ob- 

Since passively-evolving galaxies hardly accumulate stel- servational studies at higher redshift have suggested that 

lar mass only by star formation, these results suggest that more massive galaxies have older stellar population even 

some fraction of star-forming galaxies migrates into the at z ~ 1 and formed their stars earlier and more rapidly 
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than low-mass galaxies (e.g., Brinchmann & Ellis 2000; 
Juneau et al. 2005; Feulner et al. 2005; Bundy et al. 2006; 
Vcrgani et al. 2008). For quiescent galaxies with little star 
formation, Ilbert et al. (2010) and Pozzctti et al. (2010) 
reported that the number density of low-mass quiescent 
galaxies increases more rapidly than massive galaxies from 
z ~ 1 to z ~ 0. Recently, Peng et al. (2010) proposed that 
the mass-dependent quenching mechanism which ceases 
star formation preferentially for more massive galaxies is 
needed to explain no evolution of the characteristic mass 
M* for star-forming galaxies. On the other hand, the in- 
crease of low-mass quiescent galaxies at z ^ 1 is explained 
by the environmental effects in their scenario. Since the 
clustering of both bright and faint quiescent galaxies has 
been observed to be stronger than star-forming galaxies 
(e.g., Zehavi et al. 2005; McCracken et al. 2008), massive 
dark matter halos may play some role in the quenching 
of star formation. It is important to investigate how the 
evolution of quiescent galaxies and the quenching of star 
formation depend on stellar mass at higher redshift, as 
previous studies of the evolution of the global star forma- 
tion rate density and stellar mass density in the universe 
suggest that the cosmic star formation rate has a peak 
around z ~ 2 and that a significant fraction of the stellar 
mass in the present universe had been formed at z ~ 1-3 
(e.g., Hopkins & Beacom 2006; Wilkins et al. 2008). 

At z ~ 1.5-2, several studies have found massive quies- 
cent galaxies by spectroscopic observations of red objects 
(e.g., Cimatti et al. 2004; Saracco et al. 2005; Abraham 
et al. 2007; Kriek et al. 2006; Kriek et al. 2008; Kriek 
et al. 2009; Onodera et al. 2010) and by color selection 
techniques with imaging data (e.g., Daddi et al. 2004; 
Williams et al. 2009; Ilbert et al. 2010; Cameron et al. 
2010) or SED fitting techniques with multi-band pho- 
tometry (e.g., Grazian et al. 2007; Fontana et al. 2009; 
Whitaker et al. 2010). Ilbert et al. (2010) reported that 
the number density of massive quiescent galaxies with 
M star ~ 10 11 M Q rapidly increases from z ~ 1.75 and to 
z *~ 1, while it evolves only mildly at z ^ 1. Furthermore, 
Fontana et al. (2009) also found that the quiescent fraction 
in massive galaxies with M sta r > 7 x 10 10 M© significantly 
increases between z ~ 2 and z ~ 1. The low fraction of 
quiescent galaxies at z > 2 reflects that there are many 
active star-forming galaxies which dominates the massive 
galaxy population (e.g., Papovich et al. 2007; Grazian et 
al. 2007). Thus 1 < z < 2 seems to be the formation epoch 
of these massive quiescent galaxies. 

On the other hand, Kajisawa et al. (2009) found that 
the number density of ~ M* (~ 10 11 M©) galaxies evolves 
more strongly than low- mass galaxies at 1 < z J$ 3. Such 
mass-dependent evolution of the number density may in- 
dicates that the star formation history and/or stellar mass 
assembly of galaxies already strongly depended on stellar 
mass at that epoch when most stars seen in the present 
universe were formed. Therefore it is interesting to inves- 
tigate the mass-dependence of the evolution of quiescent 
galaxies up to z ~ 2. However, low-mass quiescent galaxies 
at z ^ 1-1.5 has not yet been investigated so far. Since qui- 
escent galaxies with relatively old stellar population tend 



to have high stellar mass-to- luminosity (M/L) ratios, the 
observed fluxes of quiescent galaxies with a given stellar 
mass become relatively faint (e.g., Fontana et al. 2003). 
In order to sample low-mass quiescent galaxies with high 
completeness, very deep near-infrared (NIR) data are re- 
quired (e.g., Kajisawa & Yamada 2006). 

In this paper, we study the evolution of quiescent galax- 
ies at 0.5 < z < 2.5 as a function of stellar mass using 
very deep NIR data from MOIRCS Deep Survey (MODS, 
Kajisawa et al. 2006; Ichikawa et al. 2007). The MODS 
data reach ~ 23-24 Vega magnitude (~ 25-26 AB mag- 
nitude) in the K s band, and they allow us to construct a 
stellar mass limited sample of quiescent galaxies down to 
~ 10 10 Mq even at z ~ 2. Section 2 describes the obser- 
vational data and the procedures of source detection and 
photometry. Details of the SED fitting analysis of the de- 
tected objects are given in Section 3. In Section 4, we use 
the results of the SED fitting to select quiescent galax- 
ies and construct a stellar mass-limited quiescent sample. 
We show the evolution of the number density and frac- 
tion of quiescent galaxies as a function of stellar mass in 
Section 5. In Section 6, we compare the results with pre- 
vious studies and discuss their implications. A summary 
is presented in Section 7. 

We use a cosmology with Ho=70 km s _1 Mpc -1 , O m = 
0.3 and 17a = 0.7. The Vega-referred magnitude system is 
used throughout this paper, unless stated otherwise. 

2. Observational Data and Photometry 

We use the if s -selected sample of the MODS in the 
GOODS-North region (Kajisawa et al. 2009, hereafter 
K09), which is based on our deep J_ff J^-bands imag- 
ing data taken with MOIRCS (Suzuki et al. 2008) on the 
Subaru telescope. Four MOIRCS pointings cover ~ 70% 
of the GOODS-North region (~ 103.3 arcmin 2 , hereafter 
referred as "wide" field) and the data reach J = 24.2, 
H = 23.1, K = 23.1 (5ct, Vega magnitude). One of the 
four pointings is the ultra-deep field of the MODS (~ 
28.2 arcmin 2 , hereafter "deep" field), where the data reach 
J = 25.1, = 23.7, if = 24.1. A full description of the ob- 
servations, reduction, and quality of the data is presented 
in a separate paper (Kajisawa et al. 2011). 

The source detection was performed in the if s -band im- 
age using the SExtractor image analysis package (Bertin 
& Arnouts 1996). At first, we limited the samples to 
K s < 23 and K s < 24 for the wide and deep fields, where 
the detection completeness for point sources is more than 
90% (Kajisawa et al. 2011). Then we measured the 
optical-to-MIR SEDs of the sample objects, using the pub- 
licly available multi-wavelength data in the GOODS field, 
namely KPNO/MOSAIC (U band, Capak et al. 2004), 
Hubble Space Telescope/ Advanced Camera for Surveys 
(H ST /ACS; B, V, i, z bands, version 2.0 data; M. 
Giavalisco et al. 2010, in preparation; Giavalisco et al. 
2004) and Spitzer/IRAC (3.6/xm, 4.5/im, 5.8^m, DR1 
and DR2; M. Dickinson et al. 2010, in preparation), as 
well as the MOIRCS J- and ii-bands images. Details 
of the multi-band aperture photometry are presented in 
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Kajisawa et al. (2011). Following K09, we used objects 
which are detected above 2a level in more than two other 
bands in addition to the 5a detection in the -fT^-band, be- 
cause it is difficult to estimate the photometric redshift 
and stellar mass of those detected only in one or two 
bands. The number of those excluded by this criterion 
is negligible (21/6402 and 42/3203 for the wide and deep 
fields, respectively). 

We also used public Spitzer /MIPS 24/im data in the 
GOODS-North (DR1+, M. Dickinson et al, in prepa- 
ration) to measure 24/im fluxes of the sample galaxies. 
The IRAF/DAOPHOT package (Stetson 1987) was used 
for the photometry to deal with many blended sources 
on the image properly. The DAOPHOT software fitted 
blended sources in a crowded region simultaneously with 
the PSF of the image. We used the source positions in 
the MOIRCS if s -band image as a prior for the centers 
of the fitted PSFs in the photometry For most objects, 
for which the residual of the fitting is negligible, the 5a 
limiting flux is ~ 20 //Jy. Details of the photometry on 
the MIPS image and the error estimate are described in 
Kajisawa et al. (2010) and Kajisawa et al. (2011). 

3. SED fitting analysis 

K09 estimated the redshift and stellar mass of the K s - 
selected galaxies mentioned above and constructed a stel- 
lar mass-limited sample to study the evolution of the stel- 
lar mass function. We here use the same photometric 
redshift as in K09 and carry out a similar SED fitting of 
the multi-band photometry where the grid of model pa- 
rameters is slightly changed from that in K09 in order to 
estimate age and star-formation time scale r (see below) 
optimally for the selection of quiescent galaxies. 

In K09, we show the photometric rcdshifts agree well 
with spectroscopic redshifts (Figure 1 in K09). If avail- 
able, we adopted spectroscopic redshifts from the litera- 
ture (Cohen et al. 2000; Cohen 2001; Dawson et al. 2001; 
Wirth et al. 2004; Cowie et al. 2004; Treu et al. 2005; 
Rcddy et al. 2006; Barger et al. 2008; Yoshikawa et al. 
2010), and performed the SED fitting fixing the redshift 
to each spectroscopic value for these galaxies. 

We performed the SED fitting of the multi-band pho- 
tometry (UBVizJHK, 3.6/im, 4.5/zm, and 5.8/im) with a 
population synthesis model. We adopted the GALAXEV 
population synthesis model (Bruzual & Chariot 2003) for 
direct comparison with previous studies (e.g., Grazian et 
al. 2007; Fontana et al. 2009). In the model, we assumed 
exponentially decaying star formation histories with the 
decaying timcscalc r ranging between 0.1 and 100 Gyr. 
We used Calzetti extinction law (Calzetti et al. 2000) in 
the range of E(B -V) = 0.0-1.0. Metallicity is changed 
from 1/50 to 1 solar metallicity. The model age is changed 
from 50 Myr to the age of the universe at the observed red- 
shifts. We used 18 grid of r and 60 grid of age (smaller 
number of grid at higher redshift depending on the age 
of the universe), which are roughly equally spaced on a 
logarithmic scale and are slightly finer than those used in 
K09 to select quiescent galaxies by a criterion with age/r 



(see the next section). We assume Salpeter IMF (Salpcter 
1955) with lower and upper mass limits of 0.1 and 100 M 
for easy comparison with the results in other studies. If we 
assume the Chabrier-like IMF (Chabrier 2003), the stel- 
lar mass is reduced by a factor of ~ 1.8. The redshifts 
mentioned above and stellar M/L ratios from the best- fit 
templates were used to calculate the stellar mass. The 
estimated stellar mass is reasonably consistent with that 
in K09, and the slight change in the grid of the model pa- 
rameters does not affect the estimate of the stellar mass. 
The uncertainty of the stellar mass is estimated taking 
into account of the photometric redshift error for those 
without spectroscopic redshift, and is discussed in detail 
in K09. In the next section, we will discuss the effect of 
the stellar mass error on the sample selection. 

4. Sample selection 

4-1- Selection of quiescent galaxies 

Taking advantages of the deep MOIRCS JHK s -bands 
photometry and the wide wavelength coverage of the 
multi-wavelength ancillary data of the GOODS, we used 
the results of the SED fitting described in the previous 
section to select quiescent galaxies. Several studies have 
used the SED fitting technique for the selection of qui- 
escent galaxies (Zucca et al. 2006; Arnouts et al. 2007; 
Grazian et al. 2007; Salimbcni et al. 2008; Fontana et 
al. 2009). Following Fontana et al. (2009) and Grazian 
et al. (2007), we use the criterion with age/r, assum- 
ing SFR oc exp(— age/r) in the SED fitting with the 
GALAXEV model. In this study, we adopt age/r > 6 
as in Fontana et al. (2009), which roughly corresponds 
to SFR/M s tar 10 -11 yr" 1 for relatively small r values 
(Fontana et al. 2009). It should be noted that our lower 
limit of r > 0.1 Gyr mentioned in the previous section 
leads to a lower limit of age > 0.6 Gyr for galaxies with 
age/r > 6. The upper panels of Figure 1 show the distri- 
bution of age/r as a function of stellar mass for galax- 
ies at 0.5 < z < 1.0, 1.0 < z < 1.5, and 1.5 < z < 2.5. 
Thick dashed line shows the criterion of age /r = 6 and 
thin dashed lines represent age/r = 5 and 4. The MIPS 
24/im sources (^ 20/iJy) are also shown in the figure. The 
mid-infrared detection indicates dust enshrouded star for- 
mation (or AGN activity). It is seen that at age/r ~ 4-5, 
many galaxies are detected on the MIPS 24/mi image, es- 
pecially at M st ar ^ 10 10 ' 5 M Q , while there are only a small 
number of 24/im sources at age/r > 6. The criterion of 
age/r > 5 or 4 could lead to the contamination of the 
quiescent sample with objects with some star-forming ac- 
tivities. Therefore we use the criterion of age/r > 6 in this 
paper, although we will check the results with the criteria 
of age/r > 5 and 4 in Section 6.1. For the same reason, we 
excluded the MIPS 24//m sources with age/r > 6 from the 
quiescent sample as in Fontana et al. (2009). There are 9, 
2, and 3 such 24/^m sources with age/r > 6 at 0.5 < z < 1.0, 
1.0 < z < 1.5, and 1.5 < z < 2.5, respectively. The 24/im 
fluxes of these objects could be powered by AGN (e.g., 
Nardini et al. 2008). Wc show hard X-ray sources from 
the Chandra Deep Field North catalog (Alexander et al. 
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Fig. 1. top: agc/r vs. stellar mass for ii" s -sclectcd galaxies in the MODS field in each redshift bin. The uncertainty of age/r 
is shown as error bars for galaxies with age/r > 6. Errors in age/r also include the photometric redshift error for those with no 
spectroscopic identification. Blue open squares represent objects detected in the MIPS 24/^m image with S/N > 5 (f24/jm ^ 20 
/ijy). The MIPS 24/xm sources with age/r > 6 are shown as large blue squares. Pentagons show Chandra hard-band sources from 
the CDF-North catalog (Alexander et al. 2003). Thick dashed lino shows the criterion of age/r = 6 for quiescent galaxies and thin 
dashed lines represent age/r = 5 and 4. bottom: rest-frame U — V color vs. stellar mass for the same X s -selectcd galaxies. Gray 
open squares show galaxies with age/r > 6 in the top panels. Blue larger squares and green pentagons shows the MIPS 24/^m and 
Chandra hard-band sources with agc/r > 6, respectively Horizontal dashed line represents rest U — V = 1.2. Short-dashed and 
dashed-dotted lines represent the X s -band magnitude limit for the wide (K a = 23) and deep (K a = 24) fields, respectively (see text 
for details). 



2003) in Figure 1, and only one of these 24/im sources 
with agc/r > 6 is significantly detected in the hard X-ray 
data. We also found that the other 13 24^m sources have 
the 1.6/xm bump in their near- to mid-infrared SEDs and 
do not show power-law SEDs. Therefore most of 24/im 
sources with age/r > 6 are considered to be star- forming 
galaxies, and we excluded all these 24/mi sources from 
the quiescent sample as objects with relatively active star 
formation. 

4-. 2. Limiting stellar mass for quiescent galaxies 

The iC s -band magnitude- limited sample docs not have a 
sharp limit in stellar mass even at a fixed redshift, because 
the stellar M/L ratio at the observed K s band varies with 
different stellar populations (e.g., Kajisawa & Yamada 
2006; K09). K09 investigated how the -ff s -band mag- 
nitude limit affects the stellar mass distribution of the 
sample by using the distribution of the rest-frame U — V 
color, which reflects the stellar M/L ratio well, as a func- 
tion of stellar mass (Figure 3 in K09). We here use the 
similar method to determine the limiting stellar mass for 



quiescent galaxies. The bottom panels of Figure 1 show 
the rest-frame U — V color vs. stellar mass for the K s - 
selected galaxies. Quiescent galaxies selected by age/r > 6 
in the previous section are shown as gray squares. The K s - 
band magnitude limits for the wide (K s = 23) and deep 
(K s = 24) fields are shown as Short-dashed and dashed- 
dotted lines, respectively. All objects with stellar mass 
larger than the line at a fixed U — V color (on the right 
side of the line in the figure) are brighter than the mag- 
nitude limit. We used the GALAXEV model with vari- 
ous star formation histories to calculate the line (see K09 
for details). The slope of the lines represents that galax- 
ies with redder U — V colors tend to have higher stellar 
M/L ratios and that at a fixed i4T s -band flux (luminosity), 
the limiting stellar mass for red galaxies is systematically 
higher than that for blue galaxies. 

In the figure, quiescent galaxies with age/r > 6 tend 
to have red U — V colors, typically U — V > 0.8, espe- 
cially at M st ar ^ 10 10 Mq, and form the well-known red 
sequence. It is seen that the color distribution of the qui- 
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Fig. 2. Uncertainty of the estimated stellar mass as a function of stellar mass in each redshift bin. Red circles show quiescent 
galaxies with age/r > 6 and blue dots represent those with agc/r < 6. Solid line shows the median values at each stellar mass for 
quiescent galaxies. Vertical dashed-dotted line represents the limiting stellar mass for the quiescent sample. For objects without 
spectroscopic redshift, the photometric redshift error is taken into account in the estimate of the stellar mass uncertainty. 



escent galaxies shifts to bluer with redshift. If we select 
galaxies with U- V> 0.8, U-V>0.7, and U -V > 0.6 as 
the red sequence population at 0.5 < z < 1.0, 1.0 < z < 1.5, 
and 1.5 < z < 2.5, respectively, following Whitaker et al. 
(2010) and Brammer et al. (2009), the fraction of quies- 
cent galaxies in the red sequence population is 46% (55% 
in the red sequence galaxies with M sta r > 10 11 M Q ), 36% 
(40%), and 32% (32%) for the redshift bins. Since most 
quiescent galaxies show the rest U — V ~ 1.2 or bluer, we 
adopt U — V = 1.2 (horizontal dashed line in the figure) 
as a reference to determine the limiting stellar mass for 
the galaxies in the wide and deep fields. At 1.5 < z < 2.5, 
we can nearly completely sample the galaxies down to 
- 10 10 - 5 M Q for the wide field and to - 10 10 M for the 
deep field. In the following, we use quiescent galaxies 
with M st ar > 10 10 ' 5 Mq in the wide field and those with 
M st ar > 10 10 M in the deep field as a mass-limited qui- 
escent sample at 1.5 < z < 2.5. For the 0.5 < z < 1.0 and 
1.0 < z < 1.5 redshift bins, we also limit the quiescent sam- 
ple to those with M s tar > 10 10 M© for a fair comparison 
with the 1.5 < z < 2.5 bin, while the limiting stellar mass 
is lower than 10 10 M Q in the both fields for these redshift 
bins. The number of galaxies in the mass-limited quies- 
cent sample is 130 at 0.5 < z < 1.0, 70 at 1.0 < z < 1.5 and 
50 at 1.5 <z< 2.5. 

It should be noted that our multi-band data other than 
A s -band are also deep enough to select quiescent galaxies 
down to the limiting stellar mass mentioned above. For 
example, even quiescent objects at z ~ 2-2.5 with M sta r > 
10 10 M (> 10 10 - 5 M for the wide field) whose stellar 
age is as old as the age of the universe are expected to be 
detected at least in H, 3.6/im, and 4.5/im bands at S/N 
> 2 in addition to the A s -band detection (also in J band 
for z < 2.3). Furthermore, deep H ST /ACS data are deep 
enough to detect at least quiescent objects with M s tar = 
10 10 M© and with stellar age of ~ 0.6-1.0 Gyr in V, i, and 
z bands at z < 2.5. Therefore we can judge whether the 
criteria for the quiescent population are satisfied or not for 
galaxies with M sta r ~ 10 10 M even at z ~ 2-2.5. For dusty 



galaxies, however, it could not be the case. Objects with 
large extinction are often not detected in the ACS images, 
and the uncertainty in age/r of these objects tend to be 
large. In fact, galaxies with large errorbars in the upper 
panels of Figure 1 have relatively red SEDs in the NIR- 
MIR region and are fitted with dusty model templates. 
Since the number of such objects in the quiescent sample is 
small, however, the contamination from heavily obscured 
star-forming galaxies does not seem to affect our results 
in the following. 

Figure 2 shows the uncertainty of the estimated stel- 
lar mass for the quiescent galaxies. At M sta r ^ 10 10 M Q , 
the stellar mass errors are relatively small and typically 
< 0.2 dex. Furthermore, as seen in the next section, the 
low-mass slope of the stellar mass function for quiescent 
galaxies are nearly flat or slightly positive (a ~ 0-0.6). 
Since it suggests that the number density of low-mass qui- 
escent galaxies is relatively small, the contamination from 
lower-mass galaxies could not be significant even if these 
low-mass galaxies have slightly larger stellar mass errors. 
Therefore the incompleteness and contamination due to 
the stellar mass error are expected to be negligible. 

5. Results 

The left panel of Figure 3 shows the stellar mass func- 
tion (hereafter, SMF) of quiescent galaxies selected in the 
previous section in the different redshift bins. For compar- 
ison, those for star- forming galaxies with age/r < 4 and all 
stellar mass-selected galaxies are also shown. We use the 
V max method to calculate the comoving number density. 
K09 calculated the V max for the AVband magnitude limit 
(K s = 23 for the wide field or K s = 24 for the deep field) 
by using the best-fit model template for each object (see 
Section 3.3 in K09 for details). 

It is seen that the number density of quiescent galaxies 
gradually increases from 1.5 < z < 2.5 to 0.5 < z < 1.0 
in all mass ranges. The increase of the number density 
seems to be nearly independent of stellar mass over 10 10 - 
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10 11 ' 5 M©. The number density increases by a factor of 
~ 3 from 1.0 < z < 1.5 to 0.5 < z < 1.0 and by a factor 
of ~ 10 from 1.5 < z < 2.5 to 0.5 < z < 1.0 in this mass 
range. On the other hand, the number density of star- 
forming galaxies with age /r < 4 increases only by a factor 
of ~ 1.9 from 1.0 < z < 1.5 to 0.5 < z < 1.0 and by a factor 
of - 2.8 from 1.5 < z < 2.5 to 0.5 < z < 1.0. The number 
density evolution of quiescent galaxies is clearly stronger 
than that of star-forming galaxies. 

In order to evaluate the shape of the SMF, we used the 
STY method (Sandage et al. 1979) assuming the Schechter 
function form (Schechter 1976) as in K09. We estimated 
the best-fit values of the Schechter parameters, namely, 
the characteristic mass M* and the low-mass slope a for 
both the quiescent sample and the star-forming galaxies. 
For the star-forming galaxies, we used the same limiting 
stellar mass as a function of redshift as in K09, which is 
estimated by using the 90 percentile of U — V color at 
each stellar mass for all stellar mass-selected galaxies (see 
Section 3.2 of K09 for details). Such limiting stellar mass 
is probably suitable (or conservative) for the star-forming 
galaxies, because these galaxies tend to have relatively 
low stellar M/L ratios except for extremely dusty galax- 
ies. Further details of estimating the Schechter parame- 
ters with the STY method are explained in Section 3.3 of 
K09. 

Figure 4 shows the best-fit Schechter parameters and 
their uncertainty in the M*-a plane for each redshift bin. 
The best-fit values and their uncertainty are also summa- 
rized in Table 1. For comparison, we also show the best-fit 
M* and a for all stellar mass-selected galaxies estimated 
by K09. The low- mass slope a for the quiescent sample 
is significantly flatter than those for the star-forming and 
all stellar mass-selected samples over 0.5 < z < 2.5, al- 
though the uncertainty of a for these quiescent galaxies 
is relatively large due to the small number of the sample 
and the relatively high limiting stellar mass. On the 



other hand, the characteristic mass M* for the quiescent 
sample is similar with that for the star-forming galaxies, 
while it could be slightly smaller than the star-forming 
galaxies at 1.5 < z < 2.5. The M* for the both quiescent 
and star-forming samples is lower than that for all mass- 
selected galaxies by a factor of ~ 2. This may be due 
to the effects of the combination of the (quiescent and 
star-forming) mass functions with different shapes (differ- 
ent a). K09 found a marginal upturn around ~ 10 M© 
in the stellar mass function for the stellar mass-selected 
sample, which cannot be fitted very well with the single 
Schechter function form. On the other hand, the SMF 
for the star-forming galaxies seems to have a very weak 
or no upturn, and is fitted with the Schechter function 
better than that for all mass-selected galaxies in all red- 
shift bins. The contribution of the flatter mass function 
of quiescent galaxies at M star > 10 10 M© could cause the 
upturn in the SMF for all mass-selected galaxies, which 
has been observed at z < 1 (e.g., Drory et al. 2009; Peng 
et al. 2010). Fitting the mass function with such a upturn 
with the single Schechter function may lead to the slightly 
higher M*. 

We found no significant evolution of the shape of the 
SMF for the quiescent sample. Although M* and a for 
these galaxies might become slightly lower and more pos- 
itive values respectively at 1.5 < z < 2.5, the large uncer- 
tainty prevents us from confirming this. The evolution of 
the shape of the SMF for the star-forming galaxies is sim- 
ilar with that for all mass-selected galaxies; the low-mass 
slope becomes slightly steeper with redshift at z > 1, while 
the characteristic mass does not evolve significantly. 

Figure 5 shows the fraction of quiescent galaxies at each 
stellar mass for the different redshift bins. The fraction of 
quiescent galaxies increases with stellar mass over 10 10 - 
10 1 15 M© in all redshift bins, though the uncertainty of 
the fraction at > 10 115 M© is very large due to the very 
small number of such massive quiescent galaxies in our 
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Table 1. Best-fit M* and a obtained with the STY method for quiescent and star- forming galaxies 
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sample (1-2 such galaxies in each bin). This mass depen- 
dence of the fraction of quiescent galaxies over 0.5 < z < 2.5 
is consistent with the above result that the shape of the 
SMF for these galaxies is flatter than those for the star- 
forming galaxies. On the other hand, the fraction of qui- 
escent galaxies at each mass decreases with redshift (in- 
creases with time) over 10 10 -10 n ' 5 M©. The quiescent 
fraction of galaxies with M s t a r = 10 n -10 115 Mq increases 
from - 25% at 1.5 < z < 2.5 to - 50% at 0.5 < z < 1.0, 
while that of galaxies with M star = 10 10 -10 10 - 5 M in- 
creases from < 5% to ~ 15% in the same redshift range. 
This is also consistent with the stronger evolution in the 
number density of quiescent galaxies than that of the star- 
forming galaxies mentioned above. 

6. Discussion 

6.1. Comparison with other studies 

We have investigated the evolution of the number den- 
sity and the fraction of quiescent galaxies as a function of 
stellar mass at 0.5 < z < 2.5, using the very deep NIR data 
of the MODS and the multi-wavelength ancillary data of 
the GOODS. Here we compare our results with previous 
studies in other general fields. 

Vergani ct al. (2008) investigated the evolution of early- 
type galaxies selected by spectral classification with the 
D n 4000 index up to z ~ 1.3, using the VVDS spectroscopic 
sample (Jab < 24). They found that the low- mass slope 
of the SMF for early- type galaxies is flatter (Aa ~ 1) than 
those of star-forming and overall galaxy population. The 
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Fig. 5. Fraction of quiescent galaxies as a function of stellar 
mass. Circles, pentagons, and squares represent galaxies at 
0.5 < z < 1.0, 1.0 < z < 1.5, and 1.5 < z < 2.5, respectively. 
Error bars are based on the Poisson statistics. For reference, 
the red fraction for galaxies at z ~ 0.1 from Bell ct al. (2003) 
is shown as the double-dotted dashed line. 

number density of early- type galaxies increases by a factor 
of ~ 3 between 1.0 < z < 1.3 and 0.7 < z < 1.0, while that 
of late-type galaxies increases by a factor of ~ 2 in the 
same redshift range. Our results seen in the 0.5 < z < 1.0 
and 1.0 < z < 1.5 bins are consistent with these results of 
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Number density and fraction of quiescent galaxies as a function of stellar mass for the different criteria of agc/r > 6 (left), 
5 (middle), and age/r > 4 (right). Symbols are the same as Figure 5. 



Vergani et al. (2008). 

Ilbert ct al. (2010) also studied the evolution of the SMF 
of quiescent galaxies selected by a color selection with the 
rest-frame NUV — r' and r' — J up to z ~ 2, using the 
Spitzcr/IRAC 3.6/xm-selected sample (itlab < 23.9) in the 
COSMOS 2 deg 2 field. Although their sample reach only 
to ~ 10 11 Mq at z > 1.5, they estimated the low-mass 
slope of the SMF for quiescent galaxies up to z ~ 1.2, and 
found that the low-mass slope for these galaxies gradually 
flatten with redshift and becomes a ~ at z ~ 0.8. The 
number density of these quiescent galaxies increases by a 
factor of ~ 3 from z ~ 1.2 to z ~ 0.7 and by a factor of ~ 10 
from z ~ 2 to z ~ 0.7. Such evolution of quiescent galaxies 
is consistent with our results seen in the previous section. 
We found that the flat low-mass slope for these galaxies 
continues uptoz~ 2. Ilbert et al. (2010) also pointed out 
that the evolution of the SMF for star-forming galaxies 
over 0.2 < z < 2 is relatively weak and that the number 
density of star-forming galaxies increases by a factor of 
~ 2.5 between z ~ 2 and z ~ 0.9. The star-forming galax- 
ies with age/r < 4 in our sample also show the similar 
evolution (the middle panel of Figure 3). 

Fontana ct al. (2009) investigated the evolution of the 
fraction of "red and dead" galaxies with age/r > 6 among 
massive galaxies with M star > 7 x 10 10 M Q up to z ~ 3.5, 
using the GOODS-MUSIC sample in the GOODS-South 
field. Since the survey field layout, the filter set used in 



the SED fitting, and the sample selection in this study 
are similar with those in Fontana et al. (2009), a direct 
comparison is possible. They found that the fraction of 
the "red and dead" galaxies gradually decreases with red- 
shift, and Figure 5 in Fontana et al. (2009) shows that the 
fraction is - 50% at z ~ 0.7, - 40% at z ~ 1.2, and - 20% 
at z ~ 2. If we limit our quiescent sample to those with 
M st ar > 7 x 10 10 M as in Fontana et al. (2009), the frac- 
tion of quiescent galaxies becomes 52% at 0.5 < z < 1.0, 
37% at 1.0 < z < 1.5, and 29% at 1.5 < z < 2.5. We es- 
timated the cosmic variance for these galaxies with the 
method introduced by Moster et al. (2010), using the clus- 
tering strength of rn ~ 10 Mpc for quiescent galaxies at 
1 ~ z < 2 observed by Williams et al. (2009) and Hartley 
et al. (2010). The correlation length of r ~ 10 Mpc cor- 
responds to the galaxy bias of b g ~ 2 at z ~ 0.75, b g ~ 3 
at z ~ 1.25, and b g ~ 4.5 at z ~ 2. The resulting relative 
cosmic variance for our quiescent sample is a v ~ 0.24 at 
0.5 < z < 1.0, <t v ~ 0.30 at 1.0 < z < 1.5, and a v - 0.26 at 
1.5<2<2.5. If we take account of the cosmic variance and 
Poisson errors, the quiescent fractions in this study and 
Fontana et al. (2009) agree well within the uncertainty. 

On the other hand, Grazian et al. (2007) presented the 
SMF for galaxies with age/r > 4 and those with age/r < 4 
at 1.4 < z < 2.5. In order to compare our results di- 
rectly with those in Grazian et al. (2007) and to check 
the effect of changing the selection criterion to our re- 
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suits, we estimated the SMF and the fraction of galaxies 
with age/r > 4 and 5 as a function of stellar mass, and 
compared them with those of galaxies with agc/r > 6 in 
Figure 6. Although the overall trend does not change sig- 
nificantly when the criteria of age/r > 5 and 4 are used, the 
low-mass slope becomes slightly steeper and the evolution 
of the number density becomes weaker as the selection 
threshold is loosened. This may be due to the gradual 
increase of the contamination from galaxies with star for- 
mation activities. For example, the fraction of massive 
quiescent galaxies with M sta r = 10 11 -10 115 M Q already 
reaches to ~ 50% at z ~ 2 in the case with age/r > 4. 
This is consistent with the similar number densities of 
galaxies with age/r > 4 and age/r < 4 at M sta r ~ 10 11 M© 
seen in Figure 8 of Grazian et al. (2007). 

In summary, the shape of the SMF of quiescent galax- 
ies at z <~ 1 and the number density evolution of mas- 
sive quiescent galaxies with M sta r ^ 10 11 M Q at 1 < z < 2 
seen in the previous section are consistent with the previ- 
ous studies in other general fields. Thus, the field-to-ficld 
variance does not seem to significantly affect our results. 
Furthermore, in this paper, we found that the evolution 
of the number density of quiescent galaxies with relatively 
low stellar mass (~ 10 10 TO 11 M©) is similar with that of 
massive ones, and that the flat low-mass slope of quiescent 
galaxies continues up to z ~ 2. 

6.2. Quenching of star formation and mass- dependent 
evolution of galaxies 

In Section 5, we found that the number density of qui- 
escent galaxies with age/r > 6 increases by a factor of ~ 3 
from 1.0 < z < 1.5 to 0.5 < z < 1.0 and by a factor of 
~ 10 from 1.5 < z < 2.5 to 0.5 < z < 1.0 over M star = 10 10 - 
10 11 - 5 M Q . We note that the number density of galaxies 
at 0.5 < z < 1.0 in our sample could be slightly overes- 
timated due to the known large-scale structures around 
the HDF-North (K09; Wirth et al. 2004). Nevertheless, 
the rapid increase of the number of quiescent galaxies in 
the universe seems to occur at 1 < z < 2. Furthermore, 
the evolution of the number density is nearly indepen- 
dent of stellar mass. The low-mass slope of the SMF for 
these quiescent galaxies is significantly flatter than those 
for star-forming and all stellar mass-selected galaxies over 
0.5 < z < 2.5. The fraction of quiescent galaxies for mas- 
sive galaxies is higher than that for low-mass galaxies in 
the redshift range. Here we consider that the increase 
of quiescent galaxies is caused by the cessation of star 
formation in some fraction of star-forming galaxies as in 
previous studies at z < 1 (e.g., Bell et al. 2007; Peng et 
al. 2010). Since the number density of quiescent galax- 
ies rapidly increases by a factor of ~ 3 between z ~ 2 
and z ~ 1.25 (~ 1.6 Gyr) and by a factor of ~ 3 between 
z ~ 1.25 and z ~ 0.75 (~ 2.3 Gyr), the quenching of star 
formation is expected to occur preferentially in more mas- 
sive galaxies over 1 < z < 2 in order to maintain the mass- 
dependence of the fraction of quiescent galaxies mentioned 
above. We estimated the quenching rate as a function of 
stellar mass by calculating the fraction of newly emerging 
quiescent galaxies between the redshift bins relative to the 



star-forming population including newly increased galax- 
ies at a given stellar mass range. For simplicity, other 
processes affecting this fraction such as the hierarchical 
merging is ignored. As a result, we expect that 7%, 18% 
and 29% of the star-forming galaxies with M star = 10 10 - 
10 10 - 5 M Q , 10 10 - 5 -10 n Mq and 10 u -10 115 M Q ceased 
star formation between 1.5 < z < 2.5 and 1.0 < z < 1.5, 
respectively, and 10%, 23% and 41% of the star-forming 
galaxies with the same stellar mass ranges became qui- 
escent between 1.0 < z < 1.5 and 0.5 < z < 1.0. It seems 
that more massive star-forming galaxies tend to cease star 
formation preferentially at 1 < z < 2. Recently, Peng et 
al. (2010) reported that such a mass-dependent quench- 
ing process could explain no evolution of the shape of the 
SMF for star-forming galaxies and the flatter low-mass 
slope and the similar characteristic mass for passive galax- 
ies at z < 1. Such "mass quenching" might work even at 
1 ^5 z ^ 2, when the rapid increase of the number den- 
sity of quiescent galaxies was observed. The character- 
istic mass M* of the SMF for the star-forming galaxies 
with age/r < 4 in our sample shows no significant evo- 
lution at the redshift range (middle panel of Figure 3), 
which seems to be also consistent with the prediction by 
Peng et al. (2010). 

K09 found the mass-dependent evolution of the num- 
ber density for all stellar mass-selected galaxies at z ^ 
1-1.5. The number density of galaxies with M star ~ 10 11 
Mq more strongly evolves than that for galaxies with 
Mstar ~ 10 10 Mq (right panel of Figure 3). On the 
other hand, the increase of the number density of qui- 
escent galaxies is nearly independent of stellar mass as 
discussed above. The mass dependence of the evolution 
for star-forming galaxies is also rather weak (middle panel 
of Figure 3). The mass-dependent evolution for all popu- 
lation could be explained by the differences between the 
quiescent and star-forming populations in the strength of 
the evolution of the number density and in the shape of the 
SMF. The quiescent population has the flatter low- mass 
slope of the SMF and shows the stronger evolution of the 
number density than star-forming galaxies. Therefore the 
rapid increase of quiescent galaxies contributes strongly to 
the number density evolution for galaxies with ~ 10 11 Mq, 
while their contribution is negligible for those with ~ 10 10 
Mq because of the rather low fraction of these galaxies. In 
this context, the more rapid increase of the number den- 
sity of ~ M* galaxies than lower-mass galaxies at z ^ 1- 
1.5 found in K09 might be explained by a mass-dependent 
quenching mechanism which preferentially ceases star for- 
mation in these massive galaxies. 

In Yamada et al. (2009), we found that X-ray selected 
AGNs at 2 < z < 4 are preferentially associated with 
more massive galaxies, and that ~ 30% of galaxies with 
M star > 10 11 M Q show the AGN activity with the X-ray 
luminosity larger than ~ 10 42 erg s^ 1 . Brusa et al. (2009) 
also found a similar mass-dependence of the AGN frac- 
tion at 1 < z < 4. The AGN feedback may be related with 
the mass-dependent star formation quenching. Bundy et 
al. (2008) investigated the fractions of red galaxies and 
X-ray selected AGNs in stellar mass-selected galaxies at 
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Fig. 7. Half light radius in the observed K s band vs. A" s -band magnitude for Xs-selected galaxies in each redshift bin. Gray 
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sources with age/r > 6. Solid and dashed lines represent the median values for the quiescent sample and star-forming galaxies with 
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where the PSF varies among the different pointings of MOIRCS (FWHM = 0.45-0.60 arcsec, see Kajisawa et al. 2011 for details). 
For reference, spectroscopically confirmed stars are shown as green triangles in each panel. 



0.4<z<1.4asa function of stellar mass to compare the 
star formation quenching rate with the AGN triggering 
rate. They found that the mass-dependence and normal- 
ization of the both rates agree well, although the esti- 
mated normalization of the AGN triggering rate depends 
strongly on the assumed AGN lifetime. In order to com- 
pare our results with the star formation quenching rate in 
Bundy et al. (2008), we divided the quenching rates cal- 
culated above by the time intervals between the redshift 
bins. The results are 4% Gyr" 1 , 10% Gyr" 1 , and 18% 
Gyr -1 between z ~ 1.25 and z ~ 0.75 for galaxies with 
10 10 -10 10 5 M Q , 10 10 - 5 -10 n Mq, and 10 n -10 115 M , re- 
spectively, and 4% Gyr" 1 , 11% Gyr" 1 , and 18% Gyr" 1 
between z ~ 2 and z ~ 1.25 for the same mass ranges. 
Both the mass-dependence and normalization are consis- 
tent with the results of Bundy et al. (2008), although the 
selection of quiescent galaxies is different between Bundy 
et al. (2008) (rest-frame U — B color) and this study 
(age/r from the SED fitting). Our results may also sug- 
gest that the star formation quenching rate per unit time 
at a given stellar mass seems to be roughly constant over 
1 j5 z < 2. In a future work, we will study the evolution of 
the AGN hosts in the MODS field at 1 < z < 2 and com- 
pare it with the star formation quenching rate estimated 
here. 

On the other hand, several studies reported that the 
clustering of quiescent galaxies is stronger than that of 
star-forming galaxies even at z ~ 2 (Hartley et al. 2008; 
Hartley et al. 2010; Williams et al. 2009), which implies 
that quiescent galaxies tend to inhabit more massive dark 
matter halos when their number density rapidly increases. 
Therefore the dark matter halo mass may play some role 
in the quenching of star formation. Williams ct al. (2009) 
and Hartley et al. (2010) pointed out the possibility that 
the clustering strength of quiescent galaxies at 1 < z < 2 is 
nearly independent of their rest-frame A'-band luminosity 



(approximately stellar mass), while it is observed at the 
same redshift range that the clustering strength of the 
overall galaxy population is correlated with stellar mass 
(e.g., Ichikawa et al. 2007; Foucaud et al. 2010; Wake et al. 
2010). There may be a critical halo mass above which gas 
cooling and star formation are shut down, for example, 
due to the shock heating (e.g., Dekel & Birnboim 2006; 
Cattaneo et al. 2006). In this scenario, the (stellar) mass- 
dependent quenching of star formation could be explained 
if massive star-forming galaxies tend to inhabit more mas- 
sive dark matter halos than low-mass star-forming ones, 
which is naturally predicted by theoretical models where 
the star formation activity is simply regulated by the gas 
(and dark matter) accretion rate in the halo (e.g., Bouchc 
et al. 2010). In this context, it is interesting to investi- 
gate the mass-dependence of the clustering strength for 
quiescent and star-forming samples at 1 < z < 2 more in- 
tensively. 

The quenching of star formation may be also related 
with the surface stellar mass density of galaxies. Several 
studies at low and high rcdshifts reported that the rest- 
frame color or specific star formation rate of galaxies are 
strongly correlated with their surface stellar mass density 
(e.g., Kauffmann et al. 2006; Franx et al. 2008; Williams 
et al. 2010). Most galaxies with the surface mass density 
higher than a threshold value have red color and low star 
formation activities at each redshift. The threshold sur- 
face density seems to increase with redshift gradually. At 
z ~ 2, many studies have found massive quiescent galax- 
ies with much smaller sizes (e.g., R c ~ 1 kpc) than the 
present massive ellipticals (Daddi et al. 2005; Zirm et 
al. 2007; Longhetti et al. 2007; Cimatti et al. 2008; van 
Dokkum et al. 2008; Damjanov et al. 2009; Ryan et al. 
2010). Quiescent galaxies in our sample also tend to have 
small half light radii in the observed K s band relative 
to star-forming galaxies with similar X s -band magnitudes 
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(Figure 7). Considering their higher stellar M/L ratio and 
smaller size than star-forming galaxies, these quiescent 
galaxies are expected to have higher surface stellar mass 
densities. At 1.5 < z < 2.5, there are also -fT s -bright qui- 
escent galaxies whose sizes cannot be distinguished from 
point sources. Since compact galaxies with R e ~ 1 kpc at 
z > 1 are expected not to be resolved in the MODS K s - 
band image, these galaxies may be such compact quiescent 
objects. In order to estimate the surface stellar mass den- 
sity of galaxies quantitatively and investigate the relation 
between the mass-dependent evolution of quiescent galax- 
ies and the surface stellar mass density, we need a detailed 
analysis of the surface brightness profiles which takes ac- 
count of the seeing effect or higher-resolution NIR data 
such as those with iJST/WFC3. 

7. Summary 

We have investigated the evolution of quiescent galaxies 
at 0.5 < z < 2.5 in the GOODS-North field as a function of 
stellar mass, using the very deep NIR data of the MODS 
and the public multi-wavelength data of the GOODS. We 
performed the SED fitting of the multi broad-band pho- 
tometry with the GALAXEV population synthesis model 
and selected galaxies with agc/r > 6 and without the 
MIPS 24/im detection (f^m ^ 20/xJy) as the quiescent 
population. The deep MODS data allow us to construct a 
stellar mass-limited sample of quiescent galaxies down to 
~ 10 10 M even at z ~ 2 for the first time. Main results 
in this study are as follows. 

• The number density of quiescent galaxies rapidly 
increases by a factor of ~ 3 from 1.0 < z < 1.5 to 
0.5 < z < 1.0 and by a factor of ~ 10 from 1.5 <z< 2.5 
to 0.5 < z < 1.0. On the other hand, that of star- 
forming galaxies with agc/r < 4 increases only by 
factors of ~ 2 and ~ 3 in the same redshift ranges. 

• The increase of the number density of quiescent 
galaxies is nearly independent of stellar mass over 
10 10 -10 115 M Q , and these galaxies show no signifi- 
cant evolution in the shape of the SMF. 

• The low-mass slope of the SMF for these quiescent 
galaxies is a ~ 0-0.6, which is significantly flatter 
than those of the star- forming galaxies and all stellar 
mass-selected galaxies (a~-1.3--1.5). 

• The fraction of quiescent galaxies increases with 
stellar mass over 0.5 < z < 2.5, and decreases with 
redshift over 10 10 -10 115 M Q ; the quiescent fraction 
of galaxies with KP-IO 11 - 5 M Q increases from ~ 
25% at 1.5 < z < 2.5 to - 50% at 0.5 < z < 1.0, while 
that of galaxies with 10 10 -10 10 5 Mq increases from 
< 5% to ~ 15% in the same redshift range. 

• Provided that the increase of the number density 
of quiescent galaxies is caused by the quenching of 
star formation in some star-forming galaxies, the 
quenching process is expected to be more effective 
in more massive galaxies in order to maintain the 
mass-dependence of the fraction of quiescent galax- 
ies over 0.5 < z < 2.5. 



• Since the mass dependence of the evolution of the 
number density of star-forming galaxies is also weak, 
stronger evolution of quiescent galaxies with the 
flatter shape of the SMF would cause the more 
rapid increase of the number density of galaxies 
with M star ~ 10 11 M than lower -mass galaxies at 
l<z<2. 
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